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Abstract 



In this work we consider a stochastic version of the Primitive Equa- 
tions (PEs) of the ocean and the atmosphere and establish the existence 
and uniqueness of pathwise, strong solutions. The analysis employs novel 
CO ■ techniques in contrast to previous works [18], [23] in order to handle a 

^ ' general class of nonlinear noise structures and to allow for physically rele- 

C*) , vant boundary conditions. The proof relies on Cauchy estimates, stopping 

^i ' time arguments and anisotropic estimates. 

CN 

C~-- ■ 1 Introduction 

O 

The Primitive Equations (PEs) are widely regarded as a fundamental description 
of geophysical scale fluid flows. They provide the analytical core of large General 
Circulation Models (GCMs) that arc at the forefront of numerical simulations 
K> , of the earth's ocean and atmosphere (see e.g. [47]). In view of the wide progress 

5^ ' made in computation the need has appeared to better understand and model 

some of the uncertainties which are contained in these GCMs. This is the 
so called problem of "parameterization" . Besides all of the physical forms of 
parameterization [47, 41, 40] , stochastic modeling has appeared as one of the 
major modes in the contemporary evolution of the field (see [17, 36, 37, 43, 
28, 32, 7, 50] and also [22]). In this context there is a clear need to better 
understand the numerical and analytical underpinnings of stochastic partial 
differential equations. 

In the present article we will establish the global well-posedness of the 
stochastically forced Primitive Equations of the ocean in dimension two. While 
this system has been treated in a simplified form in previous works, for the case 
of additive noise [18] and nonphysical boundary conditions [23], our aim here 
is to go further and treat a more physically realistic version of these equations 
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in the context of a multiplicative noise. In the formulation herein we face two 
new fundamental difficulties in contrast to previous work. Firstly, due to the 
imposed boundary conditions we lose higher order cancelations in the nonlin- 
ear terms. This complicates the a priori estimates which in turn prevent the 
usage of more direct compactness arguments adopted in, [8] , [23] . On the other 
hand, due to the nonlinear multiplicative noise structure, the system may not 
be transformed into a random PDE as in [18]. For this reason we are not able to 
treat the probabilistic dependence as a parameter in the problem. The analysis 
therefore requires the usage of advanced tools both from stochastic analysis, 
namely continuous time martingale theory and stopping time arguments, and 
PDE theory which we treat in detail in a separate work [21]. 

A significant literature exists concerning the Navier-Stokes equations driven 
by a multiplicative volumic white noise forcing. See [6, 48, 13, 12, 19, 34, 16, 
8, 4, 10, 33, 20]. While our point of view is similar to some of these works 
we would like to point out that the Primitive Equations, not withstanding very 
recent results on global wcll-poscdncss in 3D, arc technically more involved than 
the Navier- Stokes equations. 

This article is dedicated to Alain Bcnsoussan on the occasion of his 70th 
birthday with friendship and admiration, and, for the second author (RT), sweet 
reminiscences of many interactions, from Junior High School, to the early papers 
on stochastic partial differential equations [5], [6], on the subject of this article, 
and to many more interactions over the years. 

1.1 Presentation of the 2D Stochastic PEs 

The 2D stochastic Primitive Equations take the form 

d t u + ud x u + wd z u — vAu — fv + d x p = F u + er u (v, T)W\, (1.1a) 

d t v + ud x v + wd z v - vAv + fu = F v + <r v {v 7 T)W 2 , (Lib) 

d z p = -pg, (1.1c) 

d x u + d z w = 0, (l.ld) 

d t T + ud x T + wd z T - /xAT = F T + cr T (v, T)W 3 , (1.1c) 

P = Po(1-I3t(T-T )). (l.lf) 

This two dimensional model may be derived from the classical three dimensional 
formulation by positing invariance in one of the horizontal directions, namely 
the y— (south-north) direction. Here (v,w) — (u,v,w), T, p denote respectively 
the flow field, the temperature and the density of the fluid being modeled. The 
coefficients v, p account for the molecular viscosity and the rate of heat diffusion. 
A further parameter /, which is a function of the earth's rotation, appears in an 
antisymmetric term and is taken constant (see below). The terms F u , F v and 
Ft correspond to external sources of horizontal momentum and heat. While 
the first two terms do not usually appear in practice we retain them here for 
mathematical generality and to allow for the possible treatment, not carried out 
here, of non-homogenous boundary conditions. 
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The white noise processes Wj, the raison d'etre of the present work may be 
written in the expansions 

a u (v,T)WA 

a v (v,T)W 2 \ =a v , T (L/)^ = ^< T ([/)Fy fc . (1.2) 

,a T (v,T)W k 

The W k s may be interpreted as the time derivatives of a sequence of independent 
standard 3-D brownian motions. However, since the sample paths of brownian 
motion are nowhere differentiable we make rigorous sense of (1.1a), (1.1b) and 
(Lie) in a time integrated sense, appealing to the theory of stochastic integra- 
tion which we consider in the Ito sense. From the physical point of view these 
terms may be introduced in the model as a means to "parameterize" physical 
and numerical uncertainties. 

We consider the evolution of (1.1) over a rectangular domain Ai = (0, L) x 
(~h,0) and label the boundary I\ = (0,i) x {0}, T b = (0,L) x {-h} and 
Ti = {0,L} x (— h, 0). We posit the physically realistic boundary conditions 

<9 z v + a v v = 0, w = 0, d z T + a T T = 0, on T,, (1.3a) 

v = 0, d x T = 0, onT,, (1.3b) 

v = 0, w = 0, d z T = 0, on IV 1 (1.3c) 



The equations and boundary conditions (1.1), (1.3) are supplemented by 
initial conditions for u, v and T, that is 

u = uq, v = vq, T = To, at t = 0. (1-4) 

The Primitive equations may be derived from the compressible Navicr-Stokes 
equations with a combination of empirical observation and scale analysis. In 
particular, since deviations of the density of the fluid from a mean value are 
small at geophysical scales, the so — called Boussinesq approximation justifies 
treating the flow as incompressible. 2 Another crucial feature, that the ocean 
and atmosphere form a thin later on the earth surface leads to the hydrostatic 
approximation which reduces the third momentum equation to (1.1c). Beyond 
its obvious numerical significance, this anisotropy in the governing equations 
has many interesting theoretical consequences. We refer the interested reader 
to the classical texts [14] and [35] for an introduction from the physical point of 
view. 

Particularly in view of the numerous complications involved in extending 
the existing deterministic model to the stochastic setting we have made some 



'^Many of these boundary conditions may be non- homogenous (that is that they may include 
suitable forcing) in general. We will consider only the homogeneous case here. 

2 The Boussinesq approximation concerns the oceans. For the atmosphere we arrive at very 
similar equations by considering the pressure as the vertical coordinate, but, for the sake of 
simplicity, the emphasis here will be on the case of the oceans. 
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simplifications for the purposes of clarity of presentation. The equation (1.1) is a 
description of the earth's ocean but all of what follows can be easily extended to 
the PEs of the atmosphere or of the coupled atmosphere-ocean system (see [29]). 
We assume moreover that the /3-plane approximation is valid. This assumption, 
that the earth is locally flat, is appropriate for regional climatological studies. Of 
course, for larger scales one must include additional terms that account for the 
curvature of the earth. Since it is convenient to work in the rotating reference 
frame of the earth's surface, an additional antisymmetric term appears in the 
momentum equations. The Coriolis parameter in this term, which we denote by 
/, depends on the earth's angular velocity and the local latitude of the region 
under investigation. In the context of the /3-plane approximation, / is usually 
a linear function of y, f = /o(l + f3y). Here we take / to be constant, but once 
again the proof is easily modified to treat the more general case. 

Several other terms have been simplified or deleted which may be reintro- 
duced in their full form with no new complications to the mathematical frame- 
work or to the proof of the main theorem. We neglect the density dependence 
on the salinity of the ocean. We therefore drop the diffusion equation that ac- 
counts for variations in salt concentration in the fluid. We also ignore further, 
possibly anisotropic, diffusion terms that may appear in both the momentum 
and temperature equations to account for subgrid scale processes, the so called 
eddy diffusion terms. Finally, as noted above, we consider only the case of 
homogenous boundary conditions. 

Dating back to a series of seminal works in the early 90's [31], [30], and 
[29] a significant mathematical literature has developed around the Primitive 
Equation. In a significant breakthrough, the global well posedness in 3-D was 
established [11], [25], [26]. Subsequent work of [27] developed alternative proofs, 
which allow for the treatment of physically relevant boundary conditions. For 
the two dimensional deterministic setting we mention [38], [9] where both the 
cases of weak and strong solutions are considered. Despite these breakthroughs 
in the 3-D system, the 2-D primitive equations seem to be significantly more 
difficult mathematically than the 2-D Navier-Stokes equations. For instance, it 
is still an open problem as to whether weak solutions of the Primitive equations 
in the deterministic setting are unique. This is a classical exercise for the 2- 
D Navier Stokes equations. In any case we refer the interested reader to the 
recent survey papers [44] and [39] (appearing in [1]) which provide a systematic 
overview of deterministic theory. Note that, in regards to notational conventions 
and earlier deterministic results the present article relies heavily on this later 
work. 

While the deterministic mathematical theory is now on a firm ground the 
stochastic theory remain underdeveloped. In [23] the existence of pathwise, 
z-weak solutions was established for a simplified model with nonlinear multi- 
plicative noise and non-physical boundary conditions. A more extended system 
was considered in [18] again for the so-called z-weak solutions but with addi- 
tive noise and periodic boundary conditions. Adapting the methods of [11] the 
3-d case with additive noise and nonphysical boundary conditions was recently 
treated in [2]. 
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In contrast, beginning with the seminal work [6], extensive investigations 
for the stochastic Navier-Stokes equations have been undertaken. For weak or 
martingale solutions we mention [48], [13], [12], [19], [34] and further references 
therein. Regarding pathwise solutions we mention [16], [8], [4], [10], [33]. In 
recent joint work of the first coauthor [20] the local and global theory of pathwise 
solution in H 1 = W 1 ' 2 was established. Some of the tools and techniques 
developed in this hnal reference play a central role herein. 

In the present work we will establish the global existence and uniqueness 
of a pathwise solution to (1.1), supplemented by (1.3) for all Uq = (u,v,T) e 
(H 1 ) 3 . We conclude this introduction with an outline of the basic difficulties 
we encounter along with the main steps in the proof. 

1.2 Basic Estimates and some Difficulties Particular to 
the Stochastic Case 

The first step in the proof is to establish the local existence, up to a strictly 
positive stopping time r, of a solution U for (1.1) in U^H\ n L 2 H 2 . Here and 
throughout the rest of the work U stands for the (prognostic) unknowns in the 
problem, U = (u,v,T) = (v, T); [/(") will denote some Galerkin approximation 
of U. Having implemented a Galerkin scheme the passage to the limit is delicate 
as it is not evident a priori how to uniformly choose r > such that 

supE ( sup \U^\ 2 Hl + f |C/ (n) |^) < oo. 

n \0<t<T Jo / 

Even if such a r were to be found it would remain unclear how to infer the 
necessary sub-sequential (strong) compactness without changing the underly- 
ing stochastic basis. To overcome these difficulties we follow [20] and perform 
Cauchy type estimates for the Galerkin solutions {[/(") } n >i associated with 
(1.1) up to a carefully chosen sequence of stopping times. Since we have suffi- 
cient uniform control of the growth of U^ n ' at time zero we are able to pass to 
the limit almost surely up to a strictly positive time. Note that this stage of 
the investigation required us to establish some novel bounds on the nonlinear 
portion of the equation in H 1 (see (2.15) below) and to make careful use of the 
equivalence of some fractional order spaces. Since a significant portion of this 
analysis is non-probabilistic in character, we have separated these delicate and 
technical points to a separate work, [21]. 

With a local solution (U, t) = ((u, v, T),r) in hand, further a posteriori 
estimates are needed to preclude the possibility of a finite time blowup. In 
previous work in the deterministic setting (which corresponds to the admissible 
case, ct = 0) successive estimates on U, d z u and d x u in Lf'L 2 n L 2 !! 1 were 
conducted to finally obtain an estimate for U in LfH^. n L 2 H 2 . See [39]. For 
the present stochastic setting several difficulties emerge which prevent a trivial 
repetition of these estimates. 

The first difficulty appears when one tries to make estimates for d z u. If, on 
the one hand, we take d z of (1.1a) and then apply Ito's formula to determine 
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an evolution equation for |<9 Z «|^ 2 / M \i we encounter terms of the form 

d zzz ud z udM. 



>M 

Due to (1.3a) and (1.3c) second order terms occur on the boundary that seem 
to be intractable a priori. If, on the other hand, following [39, Section 3.3.4], 
we attempt to multiply (1.1a) by Q{—d zz u) it is not clear what the appropriate 
stochastic interpretation of du • Q(—d zz u) should be. Here Q is the orthogonal 
complement of the vertical averaging operator and is needed to get rid of the 
pressure in the governing equations (cf. (4.12) and the remarks immediately 
following) . 

To address these difficulties we introduce an auxiliary linear stochastic evo- 
lution system with a diffusion governed by the now established local solution of 
the original system. We use this system to "subtract off" the noise terms from 
(1.1a) at the cost of a number of new random terms which we must estimate. 
While we are indeed able to treat these terms, at each order our estimates re- 
quire almost sure bounds (in u) on the norms of the solution at the previous 
order. For this reason an involved stopping time argument must be employed at 
the final step. Here we make repeated use of a novel abstract result concerning 
a generic class of stochastic processes (see Proposition 5.1) which streamlines 
the analysis. 

2 Abstract Setting 

We begin with a review of the mathematical setting for the stochastic Primi- 
tive Equations and define the pathwise solutions we will consider in this work. 
The deterministic and stochastic preliminaries are treated successively. For the 
deterministic elements we largely follow [39], to which we refer the reader for 
a more detailed treatment. For more theoretical background on the general 
theory of stochastic evolution systems we mention the classical book [15] or the 
more recent treatment in [42]. 

2.1 The Hydrostatic Approximation 

The hydrostatic approximation, in concert with the incomprcssibility and the 
boundary conditions leads one to several simple observations that allow a useful 
reformulation of (1.1). This will motivate the mathematical set-up below. 

First we consider the third component of the flow w. Notice that by inte- 
grating (l.ld) and making use of the boundary condition (1.3a) for w we infer 
that 

/•O A) 

iv(x,z) = — / d z w(x,z)dz = / d x u(x,z)dz. (2-1) 



Accordingly w = w(u) is seen to be an explicit functional of u 3 . Also no- 
tice that according to the boundary conditions (1.3a), (1.3c) we impose on w, 

3 Indeed, w, p and p are called diagnostic variables in geophysical fluid mechanics. By 
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f°_ h d x udz = 0. This implies that /° A udz is constant in x and so, due to the 
lateral boundary condition (1.3b), we conclude that 

o 
udz = 0. (2.2) 

-h 

Next we consider the pressure. By integrating the hydrostatic balance equa- 
tion (1.1c) and making use of the linear dependence of the density on the tem- 
perature (l.lf) we deduce 

p s (x) -p(x,z) = d z p(x,z)dz = -gp (1 - fir(T(x,z) - T ))dz. (2.3) 



Here p s is the surface pressure, which is unknown and a function of the horizontal 
variable only. We have therefore decomposed the pressure into two components, 
the second of which couples the first momentum equation to the heat diffusion 
equation. Rearranging above and taking a partial derivative in x we arrive at 

dxP = d x p s - /3 T gpo / d x Tdz. (2.4) 

J z 

With the above considerations we now rewrite (1.1) as: 



dtu + ud x u + w(u)d z u — vAu — fv + d x p s — Prgpo / d x Tdz 



o 



(2.5a) 
=F u + <r u (v,T)W 1 , 
d t v + ud x v + w{u)d z v - i/Au + fu = F,„ + a v (v, T)W 2 , (2.5b) 

/0 pO 

d x udz, / udz = 0, (2.5c) 

d t T + ud x T + w(u)d z T - fiAT = F T + <j t {v, T)W 3 , (2.5d) 

2.2 Basic Function Spaces 

The main function spaces used are defined as follows. Take: 

H := lu=(u, v, T) e L 2 {Mf : j udz = 1 . 
We equip H with the inner product 4 

(J7,J7*):=/ v-v J AM+ f TT^dM, U = (v,T),[/» = (v*,^). 

JM JM 



opposition u, v and T are referred to as prognostic variables and are the unknowns in an 
initial value problem which we set up below. 

4 One sometimes also finds the more general definition ({7,17') := Lv ■ v'dAi + 
k f M TT"dA4 with k > fixed. This K is useful for the coherence of physical dimensions 
and for (mathematical) coercivity. Since this is not needed here we take k = 1. 
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Here and below we shall make use of the vertical averaging operator P(j> = 
r J-h 4>{z)dz and its orthogonal complement Q4> = <f> — P(j>. Note that the 
projection operator II : L 2 (A4) 3 —¥ H may be explicitly defined according to 
U n> (Qu, v, T). We also define 

V := \u= (u,v,T) EH l {M) 3 : f udz = 0, v = on T t U T b 

Here we take the inner product ((•, •)) = v{{-, -))i + M('j 0)2 where, for given 
U = (v,T),[/» = (v«,Ttt) 

((17, [/"))i := / d x v ■ d x v l + d z w ■ d z v l dM + a v f vv' dx, 
(([/,[/%:= / d x Td x T* + d z Td z T* dM + a T I TI* dx. 

J M JTi 

Note that under these definitions a Poincare type inequality \U\ < C\\U\\ holds 
for all U £ i7 1 (Al) 3 D V. Moreover the norms || • ||^i, || ■ || may be seen to be 
equivalent over all of H 1 (M) 3 . 

Even if U is very regular many of the main terms in the abstract formulation 
of (2.5) do not belong to V (see (2. 6), (2. 9), (2.10)) As such, wc shall also make 
use of some additional auxiliary spaces: 

V := | U = (u, v, T) E H\Mf : j udz = 0, v = on T t 1 , 

Z := {U = (u,v,T) e H l {Mf : v = on T;} . 

As for V we endow both spaces with the norm || • ||. One may verify that 
n : Z — > V and is continuous on H 1 (M) 3 . 

Finally we take V12) = H 2 (A4) 3 n V and equip this space with the classical 
H 2 (M) norm which we denote by | • |( 2 ). Since a considerable portion of the 
work below will consist in making estimates for the first momentum equation 
(1.1) (or equivalcntly (2.5a)) we set for simplicity 



\u\l*{m) ■= \u\, |Vu| L 2 (A1) := ||m||, \u\ H 2 {m) := |m| (2 ), 

for u G L 2 (A4) or H 1 (A4) or H 2 (M). Note that since we will always use a lower 
case u (or as needed u", u b ) for the first component of elements in the spaces 
H, V, V(2) the context will be clear. 

2.3 The deterministic framework 

The linear second order terms in the equation are captured in the Stokcs- 
type operator A which is understood as a bounded operator from V to V 
via (AU,U^) = (([/,[/")). The additional terms in the variational formulation 
of this portion of the equation capture the Robin boundary condition (1.3a). 
They may be formally derived by multiplying —uAu,—vAv,—fJ,AT in (2.5a), 
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(2.5b), (2.5d) by test functions it",u", T", integrating over M. and integrating 
by parts. We shall make use of the subspace D(A) C V(2) given by 

D(A) ={U = (v,T) eV {2 ) : <9 z v + a v v = 0, d^T + a T T = on T h 

d x T = 0onV h d z T = 0onV b }. 

On this space we may extend A to an unbounded operator by defining 

-i>QAu\ 
-uAv , U e D{A). 
-IxAT) 

Since A is self adjoint, with a compact inverse A^ 1 : H — > -D(A) we may 
apply the standard theory of compact, symmetric operators to guarantee the 
existence of an orthonormal basis {$fc}fc>o for H of cigcnfunctions of A with 
the associated eigenvalues {\k}k>o forming an unbounded, increasing sequence. 
Note that by the regularity results in [49] or [46] we have 3>fc € D(A) C VW). 
Define 

H n =.spari{$i,...,$„}. 

Take P n and Q„ = I — P n to be the projections from H onto H n and its 
orthogonal complement respectively. For m > n let P r " = P TO — P n . 

Note that in some previous works, the second component of the pressure (cf. 
(2.4) and [39, Section 2]), is included in the definition of the principal linear 
operator A. Since this breaks the symmetry of A we relegate such terms to a 
separate, lower order operator A p , which we define from V via (A P U, [/") := 
Kgpo J M J z Tdzd x u$dA4,\/U^ £ V. Taking into account the boundary conditions 
for it" on r^ (x = 0, L), this may be extended to a map A p : V — > H via 

-PrgpoQ (f z d x Tdz 


If U e D(A), A p U 6 V and we have that 

\A P U\ <c\\U\\, \\A P U\\ <c\U\ {2) . (2.7) 

We next capture the nonlinear portion of (1.1). Accordingly we define the 
diagnostic function w by setting 

w(U)=w(u)= [ d x udz, U=(u,v,T)eV. (2.8) 



For U = (v,T),C/» = (v»,T») e V we take B(U,U*) = B X {U, 0) + B 2 (U, U^) 
where 

/ Q{ud x u*)\ ( B\{u,u % y- 



Bi(t/,t7») 




B{{u,v*)\ (2.9) 
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and 

' Q(w(u)d z J)\ ( B\{u,u^ 



B 2 2 (u,J) | . (2.10) 

\Bl{u,T% 



B 2 {U,U^) := w(u)<9 z v J | 

w(u)d z T^ 

We also set -B- 7 = B{ + B 2 , j = 1,2,3. We summarize some properties of B 
needed in the sequel 

Lemma 2.1. B is well defined as a bilinear and continuous map from V x V 
to V , from V x V(o) and Vp) x V to H . Moreover B satisfies the following 
properties and estimates: 

(i) For any U,0 £V and (B{U, 11$), [/») = 0. 

(ii) ForU,U$,U b G V 

\(B{U,U*),U b )\ < c||l7||||l7*|||l7 b | 1 / 2 ||l7 l '|| 1 / 2 . ( 2 n) 

(Hi) On the other hand if we assume that U € V U* € V(2) and U G H then 
|(B(l7,^),^>|<c|[J7|[||Crl'|| 1 /a^lt|Vj»|^|. ( 2 .12) 

In particular, for U G Vt 2 ) > 

|5(£/,£/)| 2 <c||[/|| 3 |£/| (2) . (2.13) 

Also ifU= (v, T) = (u, v, T) G V(2), U^eV and U b G H, then 

\(B(U,Ut),U b )\<c\\u\\^\u\l£\\U*\\\U l '\. (2.14) 

(iv) For U G V(2), B{U) G V and satisfies the estimate 

||B([/,[/)|| 2 < C ||£/|||[/|f 2) . (2.15) 

(v) Given U, U$ G V (2) , C/ b G H 

\{B\{u,v}),u*)\ < c\u\ l l 2 \u\%)\d^\\u\ (2.16) 

\(Bl(u,J),u b }\ < clttl^llull^l^lVan^llVai^i. ( 2 .l 7 ) 

On the other hand 

\{B\{u,u*)^)\ < c|s a «||a,«»| 1 /2||a, tt i||V2| tt k| j ( 2 .is) 

\{B\{u,u*),u>)\ <c\\ut' 2 \u\%)\d^\\u\ (2.19) 
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(vi) ForU = (v,T) € D{A) 

(B 1 (u,u),-d zz u) = —r ud x u(a v u(x, 0) + d z u(x, —h)dM (2.20) 



h 



M 



which admits the estimate 
|<B 1 (u,«) ) -5„u)|<c(|«|||tt|| 2 + |a,u| 1 / a ||a,«|| 1 / a |«| 1 / 2 ||u|| 3 / 2 ). (2.21) 

The continuity properties of B as well as the basic cancellation property 
(i) are well established in the literature. The estimates (2.16), (2.17) may be 
established as for the classical Navier-Stokes systems (see, for example, [45]). 
On the other hand the estimates (2.11), (2.12), (2.14), (2.18), (2.19), (2.21) may 
be proved with anisotropic techniques. See [46] or [23]. The property (2.15), 
which is new and requires extensive computations, may be found in [21]. 

We next capture the Coriolis forcing with the bounded operator E : H — > H 
given by 

/ -Qfv \ 

EU := fu . (2.22) 

We observe that E is also continuous from V to V and that 

\EU\<c\U\, \\EU\\ < c\\U\\. (2.23) 

Finally, for brevity of notation we shall sometimes write 

N(U) = A P U + B(U, U) +EU, U e V. (2.24) 

2.4 The stochastic framework: nonlinear, multiplicative 
white noise forcing 

It finally remains to define the white noise driven terms in (1.1). To begin we 
fix a stochastic basis S := (CI, J 7 , {J-t}t>o, P, {W k }k>i), that is a filtered prob- 
ability space with {W^jfc^i a sequence of independent standard 1-D Brownian 
motions relative to the filtration J- t . In order to avoid unnecessary complica- 
tions below we may assume that Tt is complete and right continuous (see [15]). 
Fix a separable Hilbert space it with an associated orthonormal basis {efc}. We 
may formally define W by taking W = J^k W k eu- As such W is a cylindrical 
Brownian motion evolving over ii. 

We next recall some basic definitions and properties of spaces of Hilbert- 
Schmidt operators. For this purpose we suppose that X and X are any separable 
Hilbert spaces with the associated norms and inner products given by \-\x, \'\x 
and (•, -)x (•, -)xj respectively. We denote by 

L 2 (U,X) = {Re C(il,X) :J2\ Re k\x < oo}, 
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the collection of Hilbert Schmidt operators from it to X. By endowing this 
collection with the inner product 

(R, S) La (u,x) = ^(Rek, Se k ) x , 

k 

we may consider Li2(il, X) as itself being a Hilbert space. One may readily show 
that if R(U <E L 2 (il,X) and R^ G L(X,X) then indeed RWRW e L 2 {ii,X). 

Given an X-valued predictable 5 process G £ L 2 (Q; Lf oc ([0, oo),L2(il, X))) 
one may define the (Ito) stochastic integral 



M t ■- J GdW = V f G k dW k , 
Jo ,. Jo 



as a square integrable function from H, into X. Furthermore M t is an element 
of M'x, that is the space of all X-valued square integrable martingales (see 
[42, Section 2.2, 2.3]), and, as such, {M t }t>o has many desirable properties. 
Most notably the Burkholder-Davis-Gundy (BDG) inequality holds which in 
our context takes the form 



E sup 
\t'e[o,t] 



r 
GdW 

o 



x/ -° 



t \ 1/2 

L 2 (ii,X) 



< o E / \G\* X) , (2.25) 



for any t > 0, where c is here an absolute constant. 

Given any Banach spaces X and y we denote by Bnd u (X, y), the collection 
of all mappings 

* : n x [o, oo) x x ->• y, 

such that \l/ is almost surely continuous in [0, oo) x X and 

\\t>(x)\\y<c{l + \\x\\ X ), XEX, 

where the numerical constant c may be chosen independently of t and u>. If in 
addition 

\\*(x)-*(y)\\y<c\\x-y\\ x , x,y e X 

we say that "f is in Lip u (X, y). 

With these notations now in place we define 

?ct u (v,T) 

a(U)= a v (v,T) (2.26) 

cr T (v,T) 



5 For a given stochastic basis S, let <3> = [0, oo) X Q and take Q to be the cr-algcbra generated 
by sets of the form 

(s,t]x F, < s < t < oo, F G J^ s ; {0} X F, F 6 F . 

Recall that a X valued process U is called predictable (with respect to the stochastic basis S) 
if it is measurable from ($, Q) into (X, B(X)), B(X) being the family of Borel sets of X. 
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We shall assume throughout this work that 

a : n x [0,oo) x H ^ L 2 (il,H) 
such that 

If U is an i7-valued, predictable process, then 

(2.27) 
cr(U) is an L2(il, if)-valued, predictable process, 

and 

a£Lip u (H,L2{il,H))nLip u {V,L 2 {il,V))nBnd u (V,L 2 (il,D{A))). (2.28) 

Note that under the conditions imposed above the stochastic integral L a(U)dW 
may be shown to be well defined, taking values in H for any H predictable 
U 6 L 2 (rt,Lf oc ({0 7 oo);H)). Denoting Cfc(-) = er(-)efc we may interpret this 
integral in the expansion 6 

( f Qa k u (U)dW k ) 



<j(U)dW = J^ / a k (U)dW k = Y, 

k>i Jo fe>i 




t 

a k (U)dW k , 

t 
<jit(U)dW k 



1 ; <7 T {u )avv 1 

Remark 2.1. The condition (2.28) may be weakened to 

a g L J p u (i/,i 2 (il, J ff))niip„(F,L 2 (H,y))nBnd„( J D(A),L 2 (il, J D(A))) (2.29) 

m i/ie proof of local and maximal existence of solutions below (see Proposi- 
tion 3.1). However, for the proof of global existence of solutions we need the 
stronger condition (2.28). See Remark 4-2 below, for further details. Even with 
this more restrictive condition (2.28) the theory covers a physically interesting 
class of additive and nonlinear multiplicative stochastic forcing regimes relevant 
to the 'parametrization' problem discussed in the Introduction. We refer the 
interested reader to [22] for further details and examples. 

For the external forcing terms F U ,F V , Ft we let: 



F 




6 To recover the formulation of the stochastic forcings in (1.1), (1.2) we may consider the 
special case where 

o\ = when k = (mod 3) 

a v = when k = 1 (mod 3) 

cr T = when k = 2 (mod 3) 

and take Wi = £ fe W 3k e 3k , W 2 = E fc W 3k+1 e 3k+1 , W 3 = £ fc W 3k+1 e 3k+2 . 
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We assume throughout the analysis below that F is an H-valued, predictable 
process with 

FeL 2 (n;Lf oc ([0,^),H)). (2.30) 

We shall allow for the case of probabilistic dependence in the initial data Uq = 
(uo,vq,T) as well. Specifically we assume that 

U a e L 2 (ft; V) and is Jo-measurable. (2.31) 

2.5 Definition of solutions 

With the abstract mathematical definitions for each term in the original system 
now in hand we may reformulate (2.5) as an abstract evolution equation 

dU+{AU + N(U))dt = Fdt + a(U)dW, 

u(o) = u . (2 ' 32) 

More precisely we have the following basic notion of local and global pathwisc 
solutions to the above system. 

Definition 2.1 (Pathwisc Strong Solutions of the Primitive Equations). Let 
S = (f2, T , {J-t}t>o-, IP, W) be a fixed stochastic basis. Assume that F is as in 
(2.30), that Uq satisfies (2.31) and that a satisfies (2.27), (2.28). 

(i) A pair (U, r) is a local strong (pathwise) solution of (2.32) if t is a strictly 
positive stopping time and [/(■At) is a Tt adapted process in H so that 



U(-AT)eL 2 (n;C([0,oo);V)), 
C/(r)l t < T G L 2 (n;Lf oc ([0,oo);D(A))), 

and satisfies, for every t > and every U G H , 

rtAT 

(U(tAr),U)+ / (AU + N(U),U)ds 
Jo 

rt/\r ftAr 

= (U ,U}+ (F,U)ds+ (a(U),U)dW. 



(2.33) 



o Jo 



(2.34) 



(ii) Strong solutions of (2.32) are said to be (pathwise) unique up to a stopping 
time t > if given any pair of strong solutions (U 1 ,t), (U 2 ,t) which 
coincide at t = on Cl = {^(O) = U 2 (0)}, then 

P(H^(L/ 1 (tAT)-C/ 2 (tAr)) =0;Vi > 0) = 1. 

(Hi) Suppose that {r n }„>i is a strictly increasing sequence of stopping times 
converging to a (possibly infinite) stopping time £ and assume that U is a 
continuous J- t -adapted process in H . We say that the triple (U, £, {T n } ra >i) 
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is a maximal strong solution if (U, r„) is a local strong solution for each 
n and 

sup \\U\\ 2 + I \AU\ 2 ds = 00 (2.35) 

te[o,£] Jo 

almost surely on the set {£ < 00}. 

(iv) If (U, £, {r ra } n >i) is a maximal strong solution and £ = 00 o.s. t/ien toe 
say i/iai t/ie solution is global. 

We now have a complete mathematical framework and may state, in precise 
terms, the main theorem in this work: 

Theorem 2.1. Suppose that the conditions imposed in Definition 2.1 hold. 
Then there exists a unique global solution U of (2.32). 

3 Local and Maximal Existence and Uniqueness 

The proof of local and maximal existence of solutions for (2.32) makes use of 
techniques developed for the 3D Navier-Stokes Equations [20] . Since the analysis 
here is very similar on many points to [20] our treatment will be brief in some 
details. However, one crucial step, to show that the Galerkin approximations 
associated to (2.32) are Cauchy (in appropriate spaces) is quite delicate. This is 
due to stray terms that arise from the discretization which must be controlled. 
See Proposition 3.2 below. 

Proposition 3.1. Suppose that Uq, F satisfy the conditions imposed in Def- 
inition 2.1. For a we assume (2.27) and may weaken (2.28) to (2.29). Then 
there exists a unique maximal strong solution (U, £) for (2.32). Moreover, for 
any (deterministic) t > 0, 

E( sup \U\ 2 + f A \\U\\ 2 dt') <oo. (3.1) 

Y0<t'<£At Jo J 

Proof The first step in the proof, to establish certain Cauchy estimates for the 
Galerkin approximations (3.2) of (2.32) is carried out in in Lemma 3.2. For the 
details of the passage to limit we refer the reader to [20, Proposition 4.2] and 
the remarks thereafter. 

To establish local, pathwisc, uniqueness in the sense of Definition 2.1 we note 
that the estimate (2.11) of B (in dimension 2) is the same as may be achieved 
for the Navier-Stokes non-linearity in d = 3. The proof is therefore identical to 
[20, Proposition 4.1]. 

With a local strong solution in hand it remains to extend this solution to 
a maximal existence time £ as in Definition 2.1, (iii). For this point we may 
employ an argument going back to [24]. For a more recent treatment see [20, 
Lemma 4.1, 4.2, Theorem 4.1]. Since we have the cancellation property in B 
(Lemma 3.2,(i)) the bound on the weak norms up to a possible finite time blow 
up (3.1) may be established exactly as in [20, Lemma 4.2] □ 
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3.1 Local Cauchy Estimates for the Galerkin System 

We turn now to the task of estimating the difference of solutions of the Galerkin 
system associated to (2.32) at different orders. We begin by recalling some 
definitions. A ^-adapted process U^ 1 ' e L 2 (fl, C([0, oo); H n )) is a solution of 
the Galerkin system of order n for (2.32) if it satisfies: 

dU {n) + (AU {n) + P n N{U (n) ))dt = P n Fdt + P n a{U {n) )dW, 

U (n \0) = P n U o . 

Note that by the standard theory of stochastic ordinary differential equations 
one may establish the global existence of a unique solution [/'"' at each order. 
Sec e.g. [19] for details. 

Proposition 3.2. Let {EA n '} n >i be the (global) solutions of the Galerkin sys- 
tems (3.2) and suppose that there exists a deterministic constant M such that 

\\U \\ 2 < M a.s. (3.3) 

Then 

(i) there exists a stopping time T, with r > 0, a subsequence rij and a process 
U almost surely in C([0, oo); V) (~l Lf oc ([0, oo);D(A)) such that: 

lira sup \\U (n ^ - U\\ 2 + f \A(U in ^ -U)\ 2 ds = Q, (3.4) 
J^°°te[o.r] Jo 

almost surely. 

(ii) for any p > 1, there exists a sequence of Cl n . G To, with Q n . "f fi such 
that: 



supE 
J 



E n „. ( sup \\U^\\ 2 + f T \AU^\ 2 ds) 

3 V&lO,r] JO I 



< oo (3.5) 



E ( sup j|C/j| 2 + / \AU\ 2 ds ) < oo 
\*6[o,t] Jo J 



p/2 

(3.6) 



Remark 3.1. The technical condition (3.3) is needed so that we may obtain 
the uniform pathwise bound: 

supesssupf sup ||C/ (m) || 2 + f T ' n {l + \AU {m) \ 2 )ds ) < oo. (3.7) 

See (3.8) ; (3.16) below. Note however that this condition may be removed in the 
final step of the proof of the local existence. See [20, Proposition 4-2]. 
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Proof. As in previous work [20] , the proof consists in establishing the sufficient 
conditions (3.10), (3.17) for [20, Lemma 5.1] (see also related results in [33]), 
from which (i) and (ii) follow directly. The proof makes use of some delicate 
estimates present even in the deterministic case (<r = 0) that have been carried 
out in a separate work [21]. 

We assume with no loss of generality that M > 1 and consider the stopping 
times 

r,f = inf J sup \\U {n) \\ 2 + / \AU {n) \ 2 dt' > AM) . (3.8) 

*^° [t'e[o,t] Jo J 

Note that (3.8) implies that 

sup \\U {n) \\ 2 + f \AU (n) \ 2 dt' <4M, fbr0<i<r,f. (3.9) 

t'e[o,t] Jo 



We set T™ n := T n A r m ■ The first step in the proof is to perform estimates on 
jj( m ) — [/(«) w hich we denote by i?( m '") to simplify the notation below. We will 
show that 



lim su P E( sup ||i?( m '")|| 2 + / T '" lAR^^fdt) = 0, (3.10) 

n->oo m> „ \ <t'<T M Jo J 



which is the first condition required for [20, Lemma 5.1]. 

Wc fix m > n, subtract the equations for m, n, then apply A 1 ' 2 to the 
resulting system. Note that D{A 1/2 ) = V with \\U\\ 2 = \A 1/2 U\. By the Ito 
lemma we may also infer that 

d\\BS m > n) \\ 2 +2\AR {m >^\ 2 dt 

= - 2{P m N{U ( - m) ) - P n N(U {n) ),AR {m > n) )dt 

+ 2(P^F,AR<< m ^)dt (3.11) 

+ \\PmO-(U^)-P n a(U^)\\l 2iiitV) dt 

+ 2(P m a(U (m) ) - P n a{U {n) ),AR {m > n) )dW. 

We now estimate each of the terms above with a view of finally applying a 
stochastic analogue of the Gronwall inequality, [20, Lemma 5.3]. With this 
in mind fix any pair of stopping times T a ,T{, such that < t q < t\, < r^f m . 
By integrating the above system, taking a supremum over the random interval 
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[T a , T;,] and finally taking an expected values we may infer that 



E( sup ||i? (m - n) || 2 + [ " \AR (m ' n) \ 2 dt) 

<cE||fl( m ' n )(r a )|| 2 +cE f " \((P m - P n )F,AR( m ^)\dt 

+ cE f " \(P m N(U {my ) - P n N(U {n) ),AR {m ^)\dt 



+ cE / \\P m a(U {m) ) - Pno{U^)\\l 2{ixy) dt 



cE sup 

te[T a ,r b ] 



(P m a(U (m) ) - P n a{U (n) ),AR {m ^)dW 



(3.12) 



Wc begin by addressing the 'deterministic portions' of (3.12). Using the 
equivalence fractional order spaces, (2.15) and the generalized Poincare inequal- 
ity it is shown in [21], (see (3.13) in Theorem 3.1 of [21]) that: 

\(P m N{U {m) )-P n N{U {n) ),AR {m ^)\ 

<i|AR( m <">| 2 +c(l + \AU {m) \ 2 + ||L/(">|| 4 )||i?( m >">|| 2 (3 13) 

+ TTM( 1 + ll t/(n) H 2 )( 1 + l AC/( " ) | 2 )- 



We next consider the terms which arise only in the stochastic context. The Ito 
correction term may be estimated according to 



\P m a(U^)-P n a(U {n) )\\l 2 ^ v) 

<c(\\<j{U {m) )~(r{U {n) )\\ 

1 

A,, 



L 2 (il,V) 



Mu {n) )\\ 2 



L 2 (ii,V) 



<c(\\R m > n \\ 2 + -\Mu {n) )\l 2iii , H) ) 



(3.14) 



<c 



j||_ R m,n||2 + _L (l + \AU( n) \ 2 )\ 



For the second inequality we have made use of the generalized Poincare Inequal- 
ity'. The final inequality follows from (2.29). For the stochastic integral terms 



7 We use the special case ||Q„C/t|| 2 < ^-\AU S \ 2 , which holds for any t/» € D(A). 
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we apply (2.25) and deduce 

rf 



sup 

T a <t'<Tb 



(P m cr(t/ (m) ) - P n a(U M ), ABS m ' n >)dW 
<cE ( / (P m a(U^) - P n <r(UM), AR^^) 2 L2(UM) dt' 



1/2 



<cl 



||P m cr(tf<"0) - Pncr([/W)||! 2(it>y) ||i?(™^|| 2 ^' 



1/2 



<cE sup \\R^ m ^\ 
\te[T a ,T b ] 



(3.15) 



l/2> 



jT \\P m a(U^)-P n a(U^)\\ 2 L2(iiy) dt^ 

<-E( sup ||i?( m '™)|| 2 ] 

+ ce( f 6 (||i?(' m '")|| 2 + -J-(l + |^C/(")| 2 ))dt ,N ) . 

The last inequality is achieved by applying the Schwarz inequality and then 
(3.14). 

We now gather the estimates (3.13), (3.14)), (3.15) and compare with (3.12). 
Since < r a < Tb < r^f n we conclude, using (3.9) that 



E sup j|i?( m <™)|| 2 + / \AR( m ' n ^\ 2 dt 

\te[r a ,T b ] 

<cE\\R^ n ( Ta )\\ 2 



cE 



1 



(1 + \AU^\ 2 )\\BS m ^\\ 2 + 7^(1 + \AU {n) \ 2 ) + \Q n F\ 2 I dt. 

(3.16) 



Observe that the generic constant c is independent of m, n and that (3.7) (3.10) 
now follows from the stochastic Gronwall lemma. 

It remains to establish the other requirement of [20, Lemma 5.1]. In the 
present context this translates to 



limsupP sup ||t/ (n) || 2 + 

<5^0 „ \0<t'<T,"A<5 JO 



t"/\S 



\AU (n) \ 2 dt' > M] =0, (3.17) 



for every M > M. By applying Ito we infer an equation for 1 1-> ||£^™)(£)|| 2 very 
similar to (4.41), below. Since, as for the Navier-Stokes system in d = 3 (see 
(2.12)) 

\(B(U),AU)\ < ||t/|| 3/2 |AC/| 3 / 2 , U e D{A) 
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and since the A p and E terms are lower order (see (2. 7), (2. 23)), we may es- 
tablish (3.17) with a direct application of Doob's inequality exactly as in [20, 
Proposition 3.1]. With (3.17), the proof is complete. □ 

4 Global Existence 

We now implement a series of anisotropic estimates that are used to infer global 
existence. Due to the non-commutativity introduced by the physical boundary 
conditions we must first define a new variable U that satisfies a system obeying 
the rules of ordinary calculus. We are then are able to derive suitable estimates 
for u z and then u x and finally for the entire original system in V. Since the 
resulting estimates yield only pathwise (rather than moment) bounds we must 
finally recourse to some involved stopping time arguments which make essential 
use of Lemma 5.1. 

4.1 A Change of Variable to a Random PDE and some 
Auxiliary Estimates 

We consider the linear stochastic partial differential equation 

d t u - v&u + d x p s = l t < ( a u (v,T)W 1 , (4.1a) 

d t v - uAv = l t <i<r v (y,T)Wa, (4.1b) 

d t f-fiAf = l t <t:<J T (v,T)W 3 , (4.1c) 

with £ as in Proposition 3.1. This system is supplemented with the same bound- 
ary conditions as in (1.3). We posit the zero initial condition u(0) = v(0) = 
T(0) = 0. Note that the stochastic forcing terms depend on (U, £) = ((v,T),£), 
maximal strong solution solution we found for (1.1)-(1.4) in Proposition 3.1; o 
is exactly the same as appearing in (1.1) and in particular satisfies (2.28). As 
in Section 2.5, (4.1) may be formulated in an abstract form: 

dU + AUdt = l t < ( a(U)dW, 17(0) = 0. (4.2) 

We shall need the following preliminary estimates below for U. 

Lemma 4.1. There exists a unique global pathwise strong solution of (4-2) 
taking its values in D(A). Additionally for any deterministic finite time t > 0, 
we have 

E ( sup \AIJ\ 2 ] < oo. (4.3) 

\t'e[o,t] J 

Proof. We briefly outline the formal estimates that lead to (4.3). Since (4.2) is 
linear in the unknown everything, including the global existence, may be easily 
justified with a suitable Galcrkin scheme (see e.g. [19]). 



Pathwisc Solutions of the 2-D Stochastic Primitive Equations 21 



Formally then we multiply (4.2) by A and apply the Ito lemma in H to 
deduce 



d\AU\ 2 + 2\A^ 2 tJ\ 2 dt 

=21 t < e (Ar([/), AU)dW + %<z\A<x{U)\l 2{ii ,H)dt 



(4.4) 



Fixing arbitrary t > and taking a suprcmum over t' < t and then expected 
values we infer from (2.29), (4.4) and the fact that U(0) = 0, 



E sup \AU\ 2 
V'e[o,t] , 



<E sup 
t'e[o,t] 



t' hi 

(Acr(U),AU)dW 



tA£ 

n ' L 2 (il,D(A)) 1 



(4.5) 



E/ k(^)lL W .D(AlW*' 



<iE( sup \AU\ 2 )+ce[ W(U)\ 2 L2(ilMA) dt' 

1 \t'e[o,t] J Jo 

<^E( sup |yl{7| 2 )+cE/ (1 + \\Uf)dt'. 

2 y*'e[o,t] y Jo 

For the stochastic integral terms after the first inequality we apply (2.25) and 
then estimate in a similar manner to (3.15). The final inequality is a consequence 
of the assumption (2.28) imposed on a. To complete the proof we rearrange (4.5) 
and refer to (3.1) in Proposition 3.1 to conclude (4.3). □ 

We next subtract (4.1a) from (2.5) and define U = U — U. On the ran- 
dom interval [0, £) we see that U must satisfy the following partial differential 
equation (without white noise driven forcing but with random coefficients) 

^-U + All + AJU + U) + B(U + U) + E(U + U)=F. (4.6) 

dt 

Note that, in contrast to (2.32) this new system satisfies the usual rules of 
ordinary calculus. 

We may rewrite (4.6) in a form more convenient for our purposes below: 

p + AU + A p U + B(U) + EU ^ 

=F - B(U, U) - B(U, U) - B(U, U)-EU - A p U. 
By combining Lemma 4.1 with Proposition 3.1 we may directly infer that 
Lemma 4.2. For any deterministic, finite t > we have: 

E( sup \U\ 2 +f \\U\\ 2 ds) <oo (4.8) 

\0<t'<£At Jo / 
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Finally we note that the first momentum equation included in (4.7), which 
will be the focus of our attention in the subsequent sections, is given by 

dtu + ud x u + w(u)d z u — vlS.ii — fv + d x p s — flrgpo / d x Tdz 



--F u + fv + (3 T gpo \ d x fdz 

— {ud x ii + w{u)d z u) — (ud x u + w{u)d z u) — (ud x ii + w(u)d z u) 
■-F u + fv + /3 T gpo / d x fdz~(B 1 (u,u) + B 1 (u,u) + B 1 (u,u)). 



(4.9) 



Remark 4.1. We infer from (4-8) that, 



sup \U\ 2 + \\U\\ 2 ds < K\t,uj) <oo 

0<t'<£At Jo 



(4.10) 



where here and below, K,K l , denote a.s. finite constants which depend on t, 
on the data such as norms ofUo, F and on uj through these norms and though 
stochastic integral terms driven by W . 

4.2 Anisotropic Estimates 

We now turn to the estimates for d z u. 

Lemma 4.3. Let (U,£) = ((u,v,T),£) be the unique maximal strong solution 
of (2.32) guaranteed by Proposition 3.1. Then, for every t > there exists a 
finite constant K = K 2 (t,u>) < oo depending on t,ui and the data such that 

r£At 
sup |<9 2 u| 2 + / \\d z u\\ 2 ds < K 2 a.s. (4.11) 

0<f'<£At Jo 

Proof. We multiply (4.9) by —Qd zz u and integrate over the domain M.. Fol- 
lowing closely the computations in [39] we to deduce: 

2 ~dt (l 9z "l 2 + av l"li 2 (r,)) + v W d ^W 2 + va v\dxu\ 2 L 2 {Tz) 



\Pd zz u\ 2 - 


- / F u Qd zz udM 

JM 

r / r° \ 






- Prgpo 


LU, d ' {T+t)d v 


Qd zz iidM 


- f f(i 

JM 


' + v)Qd zz udM 






+ — / ud x ii[a v u(Q,x) + d z il(x 
h J M 


,-h)} 


dM 



(4.12) 



+ / {B 1 {u,u) + B 1 {u 1 u) + B 1 {u,u))d zz udM 

JM 
= J 1 + J 2 + J 3 + Ji + J5 + J6 + J7 + J 8 . 
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Here the bottom boundary is flat which causes several terms to disappear present 
in [39]. The term 

— (B 1 (ii,u),—d zz u) = / (ud x u + w{u)d z u)Qd zz udM 



largely cancels and appears as J5 due to Lemma 2.1, (vi) above. Also we observe 
that Qd x p s = which is why we multiply (4.9) by —Qd zz it rather than —d zz u. 
The first term J\ on the right hand side of (4.12) reduces to two terms at 
z = —h and that are estimated using the trace theorem: 

|Ji|<c||f7|| 2 + ^||^f. (4.13) 

The estimates for the next three terms are direct: 

|J 2 |< C |F| 2 + ^||^|| 2 , 

|J3|<c(||c>|| 2 + ||L>|| 2 ) + ^||^|| 2 

<c(||^|| 2 + |t/| 2 2) ) + ^||^u|| 2 , (4.14) 

|J 4 |< C (||C/|| 2 + ||c>|| 2 ) + ^||^|| 2 

<<\\U\\ 2 + \U\U) + ^\\dM 2 - 
For J5 we may estimate using (2.21) and Young's inequality 

|J 6 |<c(|&|||fi|| 2 + |a s fi| 1 / 2 ||M| 1 /=»|fl|i/2|| fl ||3/2 ) 

<c(|L>|||L>|| 2 + |0^| 2 / 3 |l>| 2 / 3 ||L>|| 2 ) + ^||c^|| 2 (4.15) 

< C (|L>|||L>|| 2 + |^| 2 ||L>|| 2 ) + ^||^|| 2 . 

For J 6 , (2.16), (2.19) allow 

IJel^cdfir^llfill + llfill^N^IS,^! 

<c||u|||w|( 2 )|5 zz «| (4.16) 

<c||t/|| 2 |L>| 2 2) + ^||^|| 2 . 

For J7 we estimate with (2.16) and (2.19): 

\Jr\ <c(|6| 1/2 |*|^ 2 ||fi|| + \\u\\ 1/2 \u\^\d z u\)\d zz u\ 



1(2) II "11 ~ II "ll I "1(2) 

<c\U\f 2] (\\U\\ 2 + \d z u\ 2 ) + ^-\\d z u\f 



(4.17) 
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Finally concerning Jg = (B\{u,u) + B\(u,u),d zz u) :— Jg^ + Jg,2 we estimate 

ks.il^cN^llfiir^Hfill^liil^lia^ll 

- - v (4-18) 

<c(\U\ 2 \\U\\ 2 +\\U\\ 2 \U\f 2) ) + -\\d z u\\ 2 , 
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using (2.17), and, 

IJs.al^cHtillllfillVaiiil^ll^fill 



<c||W| 2 2) + 32HM 



(4.19) 



thanks to (2.18). 

Collecting the estimates (4.13), (4.14), (4.15), (4.16), (4.17), (4.18), (4.19) 
above we may finally observe that 

-(\d z u\ 2 + a v \u\ 2 L2{ri) )+is\\d z u\\ 2 

<c(\\U\\ 2 + \U\ 2 2) )\d z u\ 2 (4-20) 

+ c(l + \U\ 2 )\\U\\ 2 + c(l + \\U\\ 2 + \m 2 )\U\ 2 i2 )+^\F\ 2 . 
We therefore conclude that 

j t (\d z u\ 2 + a v \u\ 2 L2{Ti) ) < (\d z u\ 2 + avlulli^R, + R 2 + C\F\ 2 , (4.21) 
where 

|2 , \TT\2 

(4.22) 

1(2) 



R! := \\UV + \U\ {2) 

R 2 := c(l + |L>| 2 )||L>|| 2 + c(l + ||f7|| 2 + ||l>|| 2 )|t/| 2 



and the constants c arc as in (4.20). Note that, due to (4.10) and (4.3), for all 
t > 0, there exists a constant K = K(t,u>) such that, 

/ Hjds < K (*, w) < 00 o.s. j = l,2. (4.23) 

Jo 

The (deterministic) Gronwall inequality now yields 

SUP \d z ii\ 2 < SUp (|9 z w| 2 + Qf v |M|| 2 ( r .)) 

t'e[o,T„At] t'e[o,T„At] 

< exp ( /" i?idi' j j l^uol 2 + / (i?2 + C|F| 2 )dt' j 

<K(t, W )(l + ||(7o|| 2 + / |F| 2 cit'). 



(4.24) 

Finally, returning to (4.20), integrating over [0,r„ A t], and then neglecting 
the terms |<9 2 u| 2 + a v |u|£2(ro appearing on the left hand side of the resulting 
expression, we observe that: 

/ ||cW|| 2 di'<||t/ || 2 + / (|<9 2 w| 2 i?i+i?2+c|F| 2 )di' 
Jo Jo (4-25) 

<K(t,u). 

Combining (4.24) and (4.25), completes the proof. □ 
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We next come to the estimates for d x u. Here we show 

Lemma 4.4. The hypotheses are the same as in Lemma 4-3. Then, for every 
t > 0, there exists a finite constant K = K 3 (t,u>) < oo depending on t,LU and 
the data such that 

sup \d x u\ 2 + / \\d x u\\ 2 dt' < K 3 a.s. (4.26) 

0<t'<?At Jq 

Proof. The hypotheses being the same as for Lemma 4.3, the conclusions of that 
Lemma thus hold, and in particular (4.11). 

To determine an evolution equation for \d x u\ we multiply (4.9) by —d xx u 
and integrate over M.. After some direct manipulations, this yields 

2 di'* 9 ^' 2 + ^^l 2 + ^I^IW.) 
=Pt9Po 1(1 d x {f + f)dz) O xx udM 



JM \J z J 

F u d xx udM 

(4.27) 

2f(v + v)d xx udM 

M 

+ / (B 1 (u,u) + B 1 {u,u) + B 1 (u,u) + B 1 (u,u))d xx udM 

J M 

=Ji + J 2 + J 3 + Ja + -h + Je + Ji- 

Notice that in this case the pressure term disappears by integration in z, since 
Pdxxii = 

As above the first three terms are direct 

|Ji|<c(||L>|| 2 + rii 2 ) + ^||a^|| 2 , 

|J 2 |< C |F| 2 + ^||a^|| 2 , (4.28) 

|J3l<c(||c/|| 2 + ||c/|| 2 ) + ^||a^|| 2 . 

We may handle the term J4 as in [46], however we may also directly apply 
Lemma 2.1, (2.17), (2.18) to infer 

|j 4 | ^cdiii 1 / 2 ^!! 1 / 2 ^^! 1 / 2 !!^!! 3 / 2 + la^iia^i 1 / 2 ^^!! 1 / 2 ^^!!) 

<c(H 2 ||£|| 2 |^| 2 + \d x u\ 2 \d z u\\\d z u\\) + ^\\d xU \\ 2 (4 . 29) 

<c(|c>| 2 ||L>!| 2 + ||^|| 2 )|^| 2 + ^||^|| 2 . 

The estimates (2.16) - (2.19) allow us to treat the remaining terms J5, Jq, J-j as 
well. Indeed 

IJsl^d^r^l^l^ll^ll + ll^ll 372 ^!^ 2 )!!^!! 

(4.30) 



<cj|£/|| 2 |<7| 2 2) + -Rfi|| 2 . 
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Also 



I j 6 | <cQ^ 2 \a\lg\d x <i\ + ||fi|| 1/2 |fi|$|M)IIMI 

<c|t/| 2 2) (|cX^| 2 + \d z u\ 2 ) + ^-\\d x u\\ 2 (4.31) 



]' 



<c\\U\\ 2 \U\ 2 2) + -\\d : 



,->!■ 



Finally 



\Jr\ < cdul^llwll^llull^lul^ 2 + M|||fiH 1/2 N^ a )||MI 

< 4uM(2)(\u\\\u\\ + \d x u\ 2 ) + ^\\d x u\\ 2 (4.32) 

<c\\U\\ 2 \U\l 2) + ^\\d x u\\ 2 . 
Gathering the estimates above, we conclude that: 

±\d x u\* + v\\d x u\\ 2 
dt 

< C (|L>| 2 ||L>|| 2 + ||^|| 2 )|9^| 2 (433) 

+ c(||c>|| 2 + ||l>|| 2 + ||l>|| 2 |l>| 2 2) + ||t>|| 2 |L>| 2 2) ) + c\F\ 2 

<i? 3 |a^| 2 + i?4 + c|F| 2 , 

where R 3 := c(|l>| 2 ||l>|| 2 + \\d z u\\ 2 ) and R A := c(||l>|| 2 + \\U\\ 2 + \\U\\ 2 \U\ 2 2) + 
\\U\\ 2 \U\ 2 2 ))- Dropping the term ^||9 x w|| 2 , applying the Gronwall inequality and 
then making use of the assumed bound (4.26) we infer, using (4.10), (4.3) and 
(4.11), that 



/ ,£At \ / ,fAi 

sup \d x u\ 2 <exp / R 3 dt' ) [ \d x ii \ 2 + / (R 4 + C\F\ 2 )dt' 

o<t'<£M yjo j y Jo 



I (4.34) 

<K(t,u) < +oo. 



We then integrate (4.33) from 0,£ A t and infer, using again (4.10), (4.3) and 
(4.11), that 



fA* f£At 

' BJft-fi.\ 2 + /?„ + IFI 2 W/:' 

(4.35) 

< K(t,u) < co, 



/ ||9^|| 2 dt' < ||(7 || 2 + / (i? 3 |5^| 2 + i? 4 + \F \ 2 )dt' 
Jq Jo 



where the final inequality follows from the previous bound (4.34). This com- 
pletes the proof of Lemma 4.4. □ 

Remark 4.2. With some minor modifications to the proof, Lemma 4-4 may 
established if we merely assume that, 

sup (\U\ 2 + ||L>|| 2 + \d z u\ 2 ) + / (||l>|| 2 + |L>| 2 2) + \\d z u\\ 2 )dt' < K < oo a.s. 
t'<r„ v ' h) 

(4.36) 
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On the other hand the proof of Lemma 4-3 seems to require that 

sup \U\ 2 (2) < K < oo (4.37) 

t><i 

This condition is needed to order handle Jg appearing in (4.12). The require- 
ment (4.37) is achieved due to (4.3) but at the cost of a slightly more restrictive 
condition on a, (2.28), as compared to previous work. We underline here that 
this is the only point in this work where we require the final condition in (2.28). 

Remark 4.3. We observe that the H 1 -norm \\<p\\ of a function ip is equiva- 
lent to the norm (\(p\ 2 + \d x (p\ 2 + [eJ^I 2 ) 1 ' 2 , and the H 2 - norm \f\12) of ip is 
equivalent to the norm (\\d x ip\\ 2 + \\d z <p\\ 2 + \\(p\\ 2 ) 1 ' 2 . We then infer from (4.10) 
and Lemmas 4-3 and 4-4, that u being as in these lemmas, that for every t > 0, 
there exists a constant K = K 4 (t,u>) depending on t,U) and the data, such that 

sup ||u|| 2 + / \u\ 2 2) ds < K 4 < 00 a.s. (4.38) 

0<*'<fAt JO 

4.3 Strong estimates for U 

With the above preliminaries now in hand we may now proceed to study U in 
the strong norms, the final step of the proof of global existence. 

Lemma 4.5. Suppose that < n < 00 is a deterministic constant and let 
T n < £ be the stopping time defined by 

t„ = inf \ t > : / \u\\ 2) dt' > n > A £. (4.39) 

Then, for any t > there exists a deterministic constant K = K\ (t) depending 
on n,t and the data, such that: 

e( sup \\U\\ 2 + T \AU\ 2 dA<Kl{t). (4.40) 

\0<t'<T n At Jo J 

Proof. By the Ito formula and truncation argument (see [3]) we derive an equa- 
tion for t^ \\U(t)\\: 

d\\U\\ 2 +2\AU\ 2 dt 

= (2(F - A P U - B(U) - EU,AU) + \W(U)\\l 2 ^ v) )dt (4.41) 

+ 2(A l / 2 a{U),A 1 ' 2 U)dW. 

Note that due to Proposition 3.1 this equality holds on the interval [0, £). 

Fix arbitrary stopping times < r a < 77, < T n A t. We now make estimates 
of (4.41) on this interval in order to apply the stochastic version of the Gronwall 
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lemma in [20, Lemma 5.3]. As typical, the stochastic terms are majorized by 
applying the Burkholder-Davis-Gundy inequality (2.25), 



E sup 

T a <t'<T b 



t' 

{A 1 ' 2 a{U),A l l 2 U)dW 

1 fo 

< c E (jf V^M^M^^W)^) 
<^E( sup \\U\\ 2 ) +cE / \l + \\U\\ 2 )ds. 

2 \ Ta <t'<T b J J Ta 

By applying (2.14) we may estimate the nonlinear part of the equation 
\{B(U),AU}\ < c\\uf/ 2 \u\\£\\U\\\AU\ < c\u\ 2 2) \\U\\ 2 + \\AU\ 2 

Making use of these two observations and obvious applications of Young's in- 
equality for the lower order terms (see (2.7), (2.23) ) we may estimate 



E( sup \\U\\ 2 + I " \AU\ 2 dt' 

\T a <t'<T b J Ta 



(4.42) 
< c E||C/(r a )|| 2 + c E I \l + \F\ 2 + (1 + \u\ 2 2) )\\U\\ 2 )dt' . 

The Gronwall lemma in [20] applies to real valued, non-negative processes 
X,Y,Z,R defined on an interval of time [0,T), and such that, for a stopping 
time < r < T, 

E / (RX + Z)ds < co, 
Jo 

and such that L Rds < k a.s. Assuming that, for all stopping times < r a < 
E( sup X + J " Yds) < C (E(A(r a ) + J "(RX + Z)ds 

r a <t<T b J Ta \ J Ta 

where Co is a constant independent of the choice of r a and Tb, then 

El sup 1+ / Yds) <CE(X(0)+ / Zds 
\0<t<T Jo J V Jo 

where C = C(Cq,T,K). We now just apply this lemma with r = r n ,X = 
\\U\\ 2 ,Y = \AU\ 2 ,R = c(l + \u\ 2 2) ),Z = c(l + \F\ 2 ) and the result follows. □ 

4.4 Stopping time arguments 

We now implement the stopping time arguments that, applied in combination 
with Lemmas 4.1 - 4.5, imply that £ = oo. 



Ut) 


/■tA£ 

:= sup \0 z u\ 2 + / \\d z u\\ 2 dt' 

0<t'<tA£ JO 


K*(t) 


rtA£ 

= sup \d x u\ 2 + / \\d x u\\ 2 dt' 

0<t'<tA£ Jo 


X(t): 


ntA(. 

:= sup \\U\\ 2 + / \AU\ 2 dt' 

0<t'<tA£ JO 
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We define the stochastic processes 



(4.43) 



and recall, with Lemmas 4.3 and 4.4 that X\(t) and X^if) are almost surely 
finite for all t > 0. For X(t), it follows from Lemma 4.5 that X(t) is a.s. finite 
for every t £ [0,r„] where r„ is defined by (4.39). 

We first aim to show that r„ f oo a.s. as n — > oo. Recalling that u = u + u, 
we observe that \u\ 2 2 <. < 2|u|( 2 ) + 2|w|? 2 ) and infer, with Chebyshev's inequality, 
that for any t > 0, 

P(t„ <t) <P[ / |u|( 2) ds>n J 



< P (Xi(i) + X 2 (i) > en) + -E / \u\( 2) ds. 



t 



Thanks to (4.3) this implies that 

lim P(t„ < i) < P(X x (i) + X 2 (i) = oo) = 0. 

Observing that the sequence r„ is a.s. increasing, we have 

P f lim r„ < t) = lim P(r„ < i) = 0, 

and hence r„ t °° a - s - as n — > oo. 

We now consider, for any M > 0, the stopping time 

er M =inf{r >0:X{r) > M} 

and, in view of applying Proposition 5.1 below we want to evaluate E,X(r n A 
<jm A t). To this end, we employ Lemma 4.5 and infer that 

supEX(r„ A a M At) < K*(t) < oo. 

M 

We finally conclude, by invoking Proposition 5.1, that X{t) < oo for any t > 0. 
This implies 

X(£(w)) < oo for a.a. cj e {£ < oo} (4.44) 

but since (U,£) is a maximal strong solution (cf. (2.35)), we perforce conclude 
that £ = oo a.s. The proof of Theorem 2.1 is thus complete. 
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5 Appendix I: An Abstract Stopping Time Re- 
sult 

We have made use of the following new result in the final steps of the proof 
above of global existence. 

Proposition 5.1. Fix (0, J 7 , P, {J-i}{>o), a filtered probability space. Let X : 
ft x [0, oo) — > R + U {oo} be an increasing cddldg stochastic process and define 

o-m = inf{r > : X(r) > M}. 

Suppose that there exists an increasing sequence of stopping times r„ such that 
r n 'f oo a.s. and such that for any fixed n > 0, t > 0: 

K„ t t '■= supEX(r„ A om A t) < oo. 

M 

Then, for a set Q C fl of full measure, 

X(t,uj) <oo, for allt € [0, oo), u € O. (5.1) 

Proof. It is sufficient to show that liniM->oo P(cm < i) — 0. Indeed since 

{X(t) < M} C {(Tm > t} 
and since ctjv/ is an increasing function of M, for any M' > M, 

{O-M >t}C {<T M < > t}, 

we have that 

P(X(t) < oo) = P(U M>0 {X(t) < A/}) 

< P(UAf> {cTiVf > t}) 

= lim F(o- u > t) 

M— >oo 

= lim (1-P(ct a/ <*)). 

M— >oo 

Give any M, n, observe that since X is right continuous and increasing, 

{c-m < t, r n >t} = {X(a M A t) > M, u M < t, t„ > t} 
C{X(a M At)>M,T n >t} 

C {X(a M Ar„At)> M}, 

and therefore 

P(ctm < t) < P(ctm < t, ^ > t) + P(r„ < t) 

< P(X(<r M A r„ A t) > M) + P(r n < t) 

< — E(X{o- M A T n A t)) + P(r„ < t) 
<^+P(r„<t). 
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Thus, for any fixed n and t 

lim P(ct m < t) < P(r„ < *). 

However, given the assumptions on r n , we have that 

lim P(r„ < i) = 0, 

n— K50 

which shows that X{t,ui) < oo a.s. for u £ f2. To determine the set $7 in 5.1 
and complete the proof, we observe that X is an increasing function of t and 
call, for each j £ N, Oj the set of full measure such that X(j, u) < oo, Vwe Oj. 
Then X(t, to) < oo for every t,0 < t < j, and we can take for ft, the intersection 
nj>ifij which is a set of full measure as well. □ 
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